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A THEORETICAL STUDY OF PHOTOVOLTAIC CONVERTERS 


By 

John H. Heinbockel* 

INTRODUCTION 

The research performed during the period March 1984 to January 1987 is 
summarized in the progress reports given in the Refs. [1] through [6]. The 
final phase of research was performed in the model simulation for the solar 
simulator pumped atomic iodine laser. The geometry for the steady state 
laser operation with axial lasant flow is illustrated in Fig. 1. The 
chemical kinetics for this laser are described in Refs. [7] and [8]. 

As a first approximation to the simulation of lasant flow and operation 
we have replaced the time derivatives in Ref. [7] by the material 
derivatives. 


M.)* ill + lilii (D 

Dt 3t 3z dt 

and we have assumed that the quantity — a u> represents the constant gas 

dt 

flow rate in the positive z direction. The first six of the chemical 
kinetic equations can then be written as 

3 X 3x 

— + o> — - = F..(x,t), i = 1, 2,...,6 (2) 

3t 3z 


♦Professor, Department of Mathematical Sciences, Old Dominion University, 
Norfolk, Virginia 23508. 


where x = col(x^, x 2 , x^, x^, Xg, Xg) is a column vector and 

*1 = [RI], x 2 * [R], x 3 = [R 2 ], x 4 = [I 2 ], x 5 = [I*], x g = [I]. 

_3 

Here we are using the notation [ ] to denote the concentration (cm ) of the 
chemical reactant. The synbol R denotes a radical of the perflouralkyl 
iodides (i-C^F^I, n - C^F^I, or C^I). We use the reaction rates for 
n - C^F^I and express the right hand side F^(x,t) of the system (2) as: 


F l = kjX 2 x 5 + k 2 x 2 Xg - ♦ 1 (z)x 1 - *4*2*1 


~ " k l x 2 x 5 ~ k 2 x 2 x 6 " 2k 3 x 2 ’ * C 4 X 2 X 1 


F 3 = k 3 X 2 " k 4 x l x 2 


F 4 = C 1 X 1 X 5 X 6 + C 2 X 1 X 6 + C 3 X 4 X 5 X 6 + C 4 X 4 X 6 “ *2^ x 4 


F 5 = + *2^ x 4 “ k l x 2 X 5 ' C 1 X 1 X 5 X 6 ” C 3 X 4 X 5 X 6 


■Ql x l x 5 _Q 2 x 4 x 5 ‘ r max ' ^5 “ A D X 5 


F 6 = ’M z ) x 4 + Q 1 x 1 x 5 + Q 2 X 4 X 5 + r max + Ax 5 “ C 1 X 1 X 5 X 6 


" 2c 2 X 1 X 6 " C 3 X 4 X 5 X 6 ” 2c 4 x 4 x 6 “ k 2 X 2 x 6 + k 4 x l x z V6 


where 


r max = Corp(x 5 ’ “ V 


2 


a = [2.0(10) 17 + .443x 1 ]‘ 1 
p = p + + p . 


The quantity p will be discussed in the next section. The kinetic 
coefficients and other constants are given by: 


k x = 7-9(10) ” 13 

c = 3.0(10) 11 

C 1 

= 1.0(10)" 33 

k 2 * 2.3(10) " U 

q 1 = 2.0(10) “ 16 

c 2 

= 8. 5(10) " 32 

k 3 = 2.6(10)' 12 

Q 2 = 1.9(10)' U 

c 3 

= 5.6(10)" 32 

k 4 = 3.0(10) ~ 16 

A = A D = 0 

C 4 

= 2.0(10) " 30 


Light Flux Density . We assume a one-dimensional model for the light flux 
density in which monochromatic radiation propagates along the z-axis of the 
laser. We let p + (z,t) denote the photon density propagating in the 
positive z direction and define p_(z,t) as the photon density 
propagating in the negative z direction. The differential equation for 
the photon densities are given by: 


1 9p + 
c at 


1 3p . 
c at 


+ — - op t (z,t)([l*] ■ -0]) 

3Z 2 

- — - op_(z,t)([I*] 

az 2 


(3) 


3 


where c is the speed of light in the medium. In the above equations the 

term a [1*3 — -D3 is the amplification factor for the active medium. 

2 

If Rj, denote the reflection coefficients for the faces z = 0 
and z = L, respectively, the above equations are subject to the boundary 
conditions 


P + (°.t) ■ V- ( o,t) 

P_(L,t) = R^U.t) 


(4) 


In the steady state case the solutions are denoted by p + (z) and p (z) 
and must satisfy 


p + (z)p_(z) = K q = a constant (5) 

for all values of z between 0 and L. This condition together with the 
boundary conditions requires that at z = 0, we have 

P + (0)pJO) = R x p?(0) = — P+(0)= Kq 

R 1 

and at z = L, we have 

p + (L)p_(L) = — P?(L) = R 2 p+(L)= K q 

r 2 


These conditions require that 


4 





( 6 ) 


and 



(7) 


For fixed values R^, R£ we guess at an initial value for Kg and 
integrate the system of steady state equations obtained from (2) and (3). 
This gives a calculated value of p + (L) which we compare with the 
theoretical value from (7). If these values are different we iterate on Kq 
until the final value of p+(L) agrees with the theoretical value. For 
this value of Kg and with a transmission coefficient given by t ffl = 1 - R£ 
we obtain the laser output power transmitted at the mirror where z = L. 

This power is given by 

P = e u t m cp +' L ^ ( w/cm2 ) (8) 

where is the radiation energy quantum from I*. 

The quantities ^(z) and ^(z) in the equations (2) are source 
terms related to the photodissociation rates of the chemical species. These 
functions are constructed from a function i|»(z) which is assumed to be a 
"Gaussian type curve" with a maximum value of unity at Zg^. Such a curve 
is illustrated in Fig. 2 and can be represented in the form 

F(z) = exp(-a(z-z QL ) 2 ) 


5 


The spread of this probability curve is determined by the coefficient a. 

X L 1 

It is required that when z - z n . ± we have F(z) = - (i.e. one half of 

UL o o 


its maximum value). This requires that 


1 

2 



which gives the value 


2.772 



Hence, we can write 

2 . 772 ( 1jJ0!J) 2 

' "La / 

Further, we modify this function by subtracting a constant value. The 
function. 


F(z) = exp [ - 


U-(z) = F(z) - F( 0) 

has the property that is zero at the points z = 0 and z = 2 Zq^. In order 
that the modified function have the value of unity at z * Zq L we divide by 
the appropriate scale factor and obtain 


6 



♦(I) . F(z) - F(0 > 
1 - F(0) 


We use the values 


x = 3.27 cm and z = 5.7 cm 
La OL 


and write 


* 1 (z) = C x C*^(z) 

♦ 2 (z) = € 2 C*<Kz) 

4 2 

where C* = 1.929(10) W/cm denotes the lamp concentration which is 
adjusted for the tube geometry and 5^, ^ are photodissociation rates of 
the laser gases. 

A graph of maximum power vs pressure is obtained from the numerical 
solution of the steady-state system of equations derived from (2) and (3) 
and is illustrated in Fig. 3. In Fig. 3, the curve with the circles 
represents data from the Ref. [7]. The lower three curves represent the 
numerical solution for K Q = .5(10) 25 , K Q = .5(10) 26 and K Q =.8(10) 26 . 

The computer program used to calculate the data is given in 
Appendix A. 
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Figure 3. Maximum power versus pressure. 
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APPENDIX A 
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C 
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c 

c 

c 

c 

c 

c 

601 


190 

199 


C 

C 

C 

C 
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PROGRAM LASER ( INPUT, OUTPUT, TAPE5= INPUT, TAFE6= OUTPUT, TAPES} 

MAIN PROGRAM 

CCMMCN/BLK3/B, B2, B3, C, AOO, BOO, EPSNU, OMEGA 
C0MM0N/BLK4/CHSI10, CHS I 20, AEARO, Z1BAR, LC 
C0MM0N/BLK7/ABC, COO, CO, 0MEG1, P, Rl, R2, TM 
C0MM0N/BLK8/Z0L, XLA 
REAL LC 
WRITE(6, 123} 

F0RMAT(1X, 2 OH START OF PROGRAM } 

NAMELIST/PARAM/P, 0MEG1, CON, COO, Rl, R2, LC 

CONTINUE 

READC5, FAR AM] 

I F( EOF ( 5} } 600, 601 
WR I TE ( 6 , 603} 

FORMAT (IX, 28HEND OF FILE ENCOUNTERED -STOP} 

STOP 1313 

P= PRESSURE, TORR 
0MEG1 = FLOW RATE, CM^SEC 
CCN= PEAK CONCENTRATION 

C00= INITIAL GUESS AT VALUE OF RHO-PLUS AT ZERO -- 
WHICH IS SQUARE OF (COO*R1} 

Rl = REFLECTIVITY AT LEFT END 

R2= REFLECTIVITY AT RIGHT END 

TM= TRANSMISSION COEFFICIENT =1-R2 

ZOL= POINT ALONG AXIS WHERE MAXIMUM ILLUMINATION OCCURS 
THE POINT 2* ZOL IS WHERE ILLUMINATION CUTS OFF 
THE ILLUMINATION IS BELL SHAPED ABOUT FT ZOL IN Z DIRECTION 
LC= LENGTH OF CAVITY 

XLA= WIDTH OF BELL SHAPED CURVE DEFINING LIGHT INTENSITY 
AT THE POINT OF ONE HALF MAX CONCENTRATION. 


CONTINUE 
CMIN= 1 . 0E18 
CMAX= 1 . 0E30 
TM= 1-R2 
CO=CON 
Cl 1= CON 
OMEGA =0MEG1 
XLA= 3 . 27 
ZOL= 5 . 7 
WRITE( 6, 198} 

FORMAT (///} 

WRITE(6, 199} XLA, ZOL, CON, OMEGA, COO, Rl, R2, P 

FORMAT ( IX, T5, 6HXLA = .E15.7, T30, 6HZ0L = , El 5. 7, T55, 6HC0N = , 

1 E15.7.T80, 8H0MEGA = , E15.7,/ , 

2 IX, T5, 6HC00 = ,E15.7, T30, 6H Rl = , F10. 7, T55, 6H R2 = ,F10.7, 

3 T80, ^H P = , E15 . 7 } 

SET UP COEFFICIENTS IN DIFFERENTIAL EQUATIONS 
CALL COEFFS 

CHOOSE LC SOME MULTIPLE OF .25 

OUTPUT EVERY .25 STEPS 0 .LE. Z .LE. LC 

N=i*LC 

INTEGRATE DIFFERENTIAL EQUATIONS FROM Z=0 TO Z-LC 
OUTPUT RESULTS EVERY LC'N STEPS 


702 


701 

C 

70-1 

C 

705 

C 

C 


C 

C 

100 

200 


XI = COO 

CALL INTEGCN] 

Y1 = ABC 

IFCYl.LT. 03 PER=.l 
IF(Yl.GT.O) PER= 10. 

CONTINUE 
C00= (PER3*C00 

IF(COO .LT. CM IN] STOP 5-132 
IF(COO .GT. CMAX] STOP 23-15 
X2=C00 

CALL INTEGCN] 

Y2=ABC 

IF((Y1*Y2] .LT. O] GO TO 701 

XI = COO 

Y1 = Y2 

GO TO 702 

CONTINUE 

Y1.Y2 OF OPPOSITE SIGN - USE INTERVAL HALVING TO SOLVE FOR Y 
CCO= (X1+X2]* . 5 
CALL INTEGCN] 

X3 = COO 
Y3= ABC 

IFCABSCY3] .LT. 0.001] GO TO 55 
CONTINUE 

IFCCY1*Y3] .LT. O] GO TO 705 

Y1 & Y3 ARE OF THE SAME SIGN 

XI = X3 
Y1 = Y3 
GO TO 701 
CONTINUE 

Y1 & Y3 ARE OF OPPOSITE SIGN 

v^3_ vn 

Y2 = Y3 
GO TO 701 

END 


FUNCTION CHSI1CZ] 

IMPLICIT REALa0CA-H, K, L, 0-Z] 

COMMONABLE -I^CHS 1 10, CHSI20, ABARO, Z1BAR, LC 

COMMON/BLK6/ZOL, XLA 

REAL Kl, K2, K3, K-1, LC 

IFCZ. LT. ABARO] GO TO lOO 

IFCZ . LT. Z1BAR] GO TO 200 

Z GREATER THAN Z1BAR 

CHS 11=0.0 

RETURN 

AA1 = EXFC -2 . 77*CZOL''XLA]**2] 

AA2 = EXPC -2 . 77aCCZ-Z0L]/XLA]aa2] 

FUNZ= CAA2-AA1]^C1 . -AA1] 

CHS 1 1 = CHS 1 1 Oa FUNZ 

RETURN 

END 



FUNCTION CHS I 2 ( Z } 

C IMPLICIT REAL*8 (A-H, K, L, O-Z] 

COMMON-'BLK4''CHSI10, CHS120, ABARO, Z1BAR, LC 
COMMON/ELK8/ZOL, XLA 
REAL Kl, K2, K3, K4, LC 
I F( Z . LT . ABARO] GO TO lOO 
IF(Z. LT. ZqBAR) GO TO 200 
C Z GREATER THAN Z1BAR 

100 CHS I 2 = 0.0 

RETURN 

200 AAl=EXF(-2. 77*(ZOL/XLA3**2} 

AA2=EXF(-2. 77*((Z-ZOLKXLA]**2) 

FUNZ = ( AA2 - AA1 ] ^ ( 1 . - AA1 3 
CHS 1 2 = CHS 1 20* FUNZ 
RETURN 
END 


SUBROUTINE COEFFS 

C IMPLICIT REAL* 8 (A-H, K, L, O-Z] 

C0MM0N/BLK2/K1, K2, K3, K4, Cl, C2, C3, C4, Ql, Q2, A, AD, TAUC 
C0MM0N-'BLK3-'B, B2, B3, C, AOO, BOO, EFSNU, OMEGA 
C0MM0N/BLK-1/CHSI 10, CHSI20, ABARO, Z1BAR, LC 
COMMON/BLK7/ABC, COO, CO, 0MEG1, P, Rl, R2, TM 
C0MM0N/BLK8/Z0L, XLA 
REAL Kl, K2, K3, K4, LC 

C 

C COEFFICIENTS IN THE DIFFERENTIAL EQUATIONS 

C 

OMEGA =0MEG1 
ABARO=0. 0 

C ABARO = START OF ILLUMINATION 

C Z1BAR=2*Z0L = POINT ON AXIS WHERE ILLUMINATION STOPS 

CHSI10=(3.04E-3}*C0 
CHS 1 20= ( 3 . 38E-2]*CO 
Z1BAR=2*Z0L 
EPSNU=1 . 5E-19 

C WATTS--' CM* CM 

A00= 2 . 0E17 
B00= . -H3 
Kl = 7 . 9E-13 
K2 = 2 . 3E-11 
K3= 2 . 6E-12 
K4 = 3 . OE-16 
C = 3 . 0E10 
Ql=2. OE-16 
Q2= 1 . 9E-11 
B=P*(3. 5E16] 

Cl=l . OE-33 
C2= 8 . 5E-32 
C3= 5 . 6E-32 
C-1 = 2 . 0E-30 
A= 0 

AD= 1 . 2E- 3 
B2=B*B 
B3 = B2*B 
RETURN 
END 


non 


SUBROUTINE FUN(2, Y, F] 

C IMPLICIT REAL* 8 (A-H, K, L, 0-2] 

DIMENSION Y(7] , F(7] 

C0MM0N/BLK1/X7, POWER 
EXTERNAL CHS 1 1 , CHS 1 2 

COMMONABLE PAKl, K2, K3, K4, Cl, C2, C3, C4, Ql, Q2, A, AD, TAUC 
C0MM0N/BLK3/B, B2, B3, C, AOO, BOO, EPSNU, OMEGA 
C0MM0N/BLK4/CHSI10, CHS 1 20, ABARO, 21 BAR, LC 
COMMONABLE 7/ ABC, COO, CO, 0MEG1, P. Rl, R2, TM 
REAL Kl, K2, K3, K-1, LC 
C 

X8=COO/(Y(7]*B] 

S1G= 1 . a ( AOO+BOO*B*Y ( 1] ] 

X7STAR= Y( 7] *B+X8 
DIF=Y(5]-. 5*Y( 6] 

F( 1] = K1*B*Y( 2]*Y( 5] -»K2*B*Y( 2]*Y( 6] -CHS 1 1 ( 2]*Y( 1] -K-1*B*Y (1]*Y(2] 

F( 2] =CHS1 1 (2]*Y( 1] -K1*B*Y( 2]*Y( 5] -K2*B*Y ( 2] *Y( 6] -2*K3*B*Y( 2] *Y( 2] 
1 -K-1*B*Y(1]*Y(2] 

F( 3] = K3*B*Y( 2]*Y( 2] -»K4*B*Y ( 1]*Y( 2] 

A1 = C1*B2*Y ( 1 ] *Y( 5] *Y (6] + C2*B2*Y ( 1 ] *Y( 6] *Y (6] *C3*B2*YM]*Y( 5] *Y( 6] 
A2=C4*B2*Y(4]*Y(6]*Y(6]-CHSI2(2]*Y(-1] 

F(-1] = A1+A2 

A3 = CHS I1(2]*Y(1] +CHSI 2( 2] *Y(-1] -K1*B*Y (2]*Y(5] 

A-1 = -C1*B2*Y(1]*Y(5]*Y(6]-C3*B2*Y(-1]*Y(5]*Y(6]-Q1*B*Y(1]*Y(5] 

A5 = -Q2*B*YM]*Y(5]-C*SIG*X7STAR*DIF-A*Y(5] ~AD*Y(5] 

F( 5] = A3+A4+A5 

A6=CHSI2(2]*Y(4] +Q1*B*Y( 1 ] *Y( 5] +Q2*B*Y( 4] *Y( 5] 
A7=C*SIG*X7STAR*DIF+A*Y(5] -C1*B2*Y( 1 ] *Y( 5] *Y( 6] 
A8=-2*C2*B2*Y(1]*Y(6]*Y(6]-C3*B2*Y(4]*Y(5]*Y(6] -AD*Y(6] 
A9=-2*C4*B2*Y(4]*Y(6]*Y(6]-K2*B*Y(2]*Y(6] * K4*B*Y( 1 ] *Y( 2] 

F( 6] = A6+A7+A8+A9 
DO lO 1 = 1, 6 

10 F(I]=F(I] AOMEGA 

F( 7] =Y(7]*D1F*B*SIG 

RETURN 

END 


SUBROUTINE INTEG(N] 

C IMPLICIT REAL*8(A-H, K, L, 0-2] 

DIMENSION Y( 7] , F( 7] , Y0( 7] , X( 7] , AA1 ( 7] , AA2( 7] , AA3( 7] , AA-1 ( 7] , U( 7] 
COMMONABLK1AX7, POWER 

COMMONABLK3AB, B2, B3, C, AOO, BOO, EPSNU, OMEGA 
COMMONABLK^acHSIIO, CHS I 20, ABARO, 21 BAR, LC 
C0MM0NABLK7AABC, COO, CO, 0MEG1 , P, Rl, R2, TM 
REAL Kl, K2, K3, K-1, LC 

INTEGRATE SYSTEM FROM 2=0 TO 2=LC USING RUNGE-KUTTA METHOD 

55 CONTINUE 

STEP = LC A FLOAT ( N ] 

NTIME= 500 

H= STEPAFLOAT ( NT I ME] 

NPR1NT= 500 

C INITIAL CONDITIONS 

20 = 0.0 
Y0( 1] =1.0 
DO 9 1 = 2, 6 
9 Y0( I ] = O . O 


c 

191 

300 

C 

10 

199 

C 

C 

100 

11 

12 

13 

11 

15 

16 

17 

18 


GUESS AT INITIAL CONDITIONS FOB X(7) AND X(8) 

X70=SQRT(Rl*COO) 

YO( 7} = X?0< / 'B 
WRITE(6, 191) 

FORMAT T7, 1HZ, T20, 4HX(1), T32, 1HX(2) , T15, 1HX(3), T57, 1HX(i) , 
1 T69, 1HX(5), T80, 1HX(6), T91, 1HX(7), T103, 4HX(8) ) 

CONTINUE 
11 = 0 

CALL FUN(ZO, YO, F) 

PRINT OUT 
DO lO 1 = 1, 7 
X(1)=B*Y0(I) 

X8=COO/X(7) 

WRITE(6, 199) ZO, (X(I), 1 = 1, 7),X8 
F0RMAT(1X, El 2. 5, 8E12. 5, El 2. 5 ) 

INTEGRATE USING STEP SIZE H XTIMES THEN PRINT OUT AGAIN 
11=11+1 

CALL FUN(ZO, YO, F) 

DO 11 1 = 1, 7 
AA1(I)=H*F(I) 

DO 12 1 = 1, 7 

U( I ) = Y0( I) + . 5*AA1 ( I ) 

Z1=Z1+.5*H 
CALL FUN(Z1, U, F) 

DO 13 1 = 1,7 
AA2( I ) =H*F( I ) 

DO 11 1 = 1, 7 

U( I ) = Y0( I ) + . 5*AA2( I ) 

CALL FUN(Z1, U, F) 

DO 15 1 = 1, 7 
AA3( I ) =H*F( I ) 

Zl=ZO+H 

DO 16 1 = 1,7 

U( I ) =Y0( I ) +AA3( I ) 

CALL FUN(Z1,U, F) 

DO 17 1 = 1,7 
AA1(I)=H*F(I) 

DO 18 1 = 1,7 

Y( I ) = Y0( I ) + ( AA1 ( I ) +AA1 ( I ) + 2* ( AA2 ( I ) + AA3( I ) ) )'6 . 

Z=ZO+H 



C UPDATE ZO, YO VALUES 

ZO=Z 

DO IS 1=1,7 
19 YO(I]=Y(I] 

IF( 1 1 . GE.NFR1NT]G0 TO 200 
GO TO 100 
200 CONTINUE 

X( 73 = B*YO( 73 
X8=COO/X(7] 

IF( (Z0-» . 5*H] . GE. LC] GO TO 500 
GO TO 300 
500 CONTINUE 

CALL FUN(ZO, YO, F] 

DO 110 1 = 1, 7 
110 X(I]=B*YO(I] 

XX7L=B*Y(7] 

X8=COO/X(7] 

RH0FL=X70/SQRT(R1*R2] 

C RHOPL= RHO-PLUS AT Z = L THEORETICAL VALUE 

C XX7L= CALCULATED VALUE OF RHO-PLUS AT Z= L 

C ABC = D I FFERENCE= XX7L-RH0PL 

DIF= ( ( (XX7L-RHOFL]/RHOPL]*10C] 

ABC=DIF 

WRITE( 6, 202]DIF, RHOPL, XX7L, COO 

202 FORMAT( IX, 1 3HDI FFERENCE = , E18.S, 2X, 1 2HTHE0RETI CAL= , E18.S, 2X, 

1 10H ACTUAL = , El 8. 9, 2X, 6HC00 = , E18.8 ] 

CALL OUTFUT(YO( 7} , ZO] 

RETURN 

END 


SUBROUTINE OUTPUT (YY, ZZ] 

COMMON/BLK3/B, B2, B3, C, AOO, BOO, EPSNU, OMEGA 
COMMON/BLK^CHSIIO, CHS120, ABARO, Z1BAR, LC 
COMMON/BLK7/ABC, COO, CO, 0MEG1 , P, Rl, R2, TM 
REAL LC 
XX7L=B*YY 

PGWER= EPSNU*TM*C*XX7L 
WRITE(6, 193]R1, R2, POWER, TM, ZZ 

193 F0RMAT(1X, 5HR1 = , F10 . 7, IX, 5HR2 = , F10 . 7, IX, 7HP0WER = .E18.10, 
1 IX, 5HTM = , F10. 8, IX, 4HL = .F15.7 ] 

RETURN 

END 


